REGULARITY "IN LARGE" FOR THE 'SD SALMON'S PLANETARY 
GEOSTROPHIC MODEL OF OCEAN DYNAMICS 
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Abstract. It is well known, by now, that the three-dimensional non-viscous planetary geostrophic model, 
with vertical hydrostatic balance and horizontal Rayleigh friction, coupled to the heat diffusion and 
transport, is mathematically ill-posed. This is because the no-normal flow physical boundary condition 
implicitly produces an additional boundary condition for the temperature at the literal boundary. This 
additional boundary condition is different, because of the Coriolis forcing term, than the no heat flux 
physical boundary condition. Consequently, the second order parabolic heat equation is over determined 
with two different boundary conditions. In a previous work we proposed one remedy to this problem 
by introducing a fourth-order artificial hyper-diffusion to the heat transport equation and proved global 
regularity for the proposed model. Another remedy for this problem was suggested by R. Salmon by 
introducing an additional Rayleigh-like friction term for the vertical component of the velocity in the 
hydrostatic balance equation. In this paper we prove the global, for all time and all initial data, well- 
posedness of strong solutions to the three-dimensional Salmon's planetary geostrophic model of ocean 
dynamics. That is, we show global existence, uniqueness and continuous dependence of the strong solutions 
on initial data for this model. 
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1. Introduction 

The starting point in the derivation of the ocean circulation models is Boussinesq equations which are 
the Navier-Stokes equations with rotation and a heat transport equation. The global existence of strong 
solution to the Navier-Stokes equations, which is a particular case of the Boussinesq equations when the 
temperature is identically zero, is one of the most challenging problems in applied analysis. However, 
geophysicists take advantage of the shallowness of the oceans and the atmosphere and introduce the 
hydrostatic balance approximation in the vertical motion. This in turn simplifies the Boussinesq model, 
and leads to the primitive equations of ocean and atmosphere dynamics (see, e.g., [14], [15], [16], [18], 
[20], [22], [34] and references therein). Further, horizontally, approximations based on the fast rotation of 
the earth, and the shallowness of the atmosphere and ocean imply the smallness of the Rossby number, 
which consequently lead to the geostrophic balance between the Coriolis force and the horizontal pressure 
gradient (cf. e.g., [11], [18], [22], [34] and references therein). By taking advantage of these assumptions and 
other geophysical considerations several intermediate models have been developed and used in numerical 
studies of weather prediction, long-time climate dynamics and large scale ocean circulation dynamics (see, 
e.g., [2], [3], [6], [7], [18], [20], [23], [26], [27], [28], [29], [36] and references therein). 

The planetary geostrophic (PG) model, the inviscid and adiabatic form of "thermocline" equations, of 
large scale ocean circulation are derived by standard scaling analysis for gyre-scale oceanic motion (see 
[17], [19] , [21], [22], [34] and [35]). They are given in their simplest dimensionless /3— plane mid-latitude 



Date: December 27, 2010. 

1 



2 



C. CAO AND E.S. TITI 



approximation by the system of equations: 

P.~fv^O, py + fu = 0, p,-T = 0, (1) 

Ux + Vy + Wz = (2) 

dtT + uTx + vTy + wT, = KyTzz , (3) 



in the domain = {{x,y,z) : {x,y) G M C M^, and z G (— /i, 0)}. Here {u,v,w) denotes the velocity 
field, p is the pressure, and T is the temperature, which are the unknowns. / = /o + /3y is the /3— plane 
mid-latitude approximation of the Coriolis force. The first two equations in (1) represent the geostrophic 
balance and the third equation represents the hydrostatic balance. The diffusive term, KyT^z is a leading 
order approximation to the effect of macro-scale turbulent mixing. Based on physical ground Samelson 
and Vallis [26] have argued that in closed ocean basin, with the no-normal-flow boundary conditions, 
this model can be solved only in restricted domains which are bounded away from the lateral boundary, 
dM X (— /i, 0). Thus, it cannot be utilized in the study of the large-scale circulation. Furthermore, it 
has been pointed out numerically in [8] that arbitrarily small linear disturbances (disturbances that are 
supported at small spatial scales) will grow arbitrarily fast when the flow becomes baroclinically unstable. 
This nonphysical growth at small scales is a signature of mathematical ill-posedness of this model near 
unstable baroclinic mode. Therefore, Samelson and Vallis proposed in [26] various dissipative schemes 
to overcome these physical and numerical difficulties. In particular, they propose to add either a linear 
Rayleigh-like drag/friction/damping or a conventional eddy viscosity to the horizontal components of the 
momentum equations, and a horizontal diffiision in the thermodynamic equation (subject to no-heat-flux 
at the lateral boundary.) The planetary geostrophic (PG) model with conventional eddy viscosity has 
been studied mathematically in [4], [24], [25]. In [4] we show the global existence and uniqueness of weak 
and strong solutions to this 3D viscous PG model. We also provide rigorous estimates, depending on the 
various physical parameters, for the dimension of its global attractor. In the case where the dissipative 
scheme for the horizontal momentum is the linear drag Rayleigh friction it is observed that the second 
order parabolic PDE that governs the temperature (the thermodynamic equation) has, due to the Coriolis 
force, too many boundary conditions to be satisfied, and hence it is over determined and is ill-posed (see, 
e.g., the discussion in section 2 below, [5], [26] and the references therein). To remedy this situation it is 
argued in [26] that one would have to add to the thermodynamic equation a higher order (biharmonic) 
horizontal diffusion in order to be able to satisfy both physical boundary - the no-normal-flow and no- 
heat-flux boundary conditions - (cf. e.g., [5], [26], [27]). In [5], we introduce, instead, a new PG model 
with an appropriate artificial horizontal "hyperdiffusion" term, to the heat equation, which involves the 
Coriolis parameters. Under the two natural physical boundary conditions we are able to prove in [5] the 
global existence and uniqueness of the strong solutions. Moreover, we also show the existence of the flnite 
dimensional global attractor. 

To overcome the above mentioned non-physical baroclinical instabilities and numerical ill-posedness 
Salmon introduced in [22] the following alternative planetary geostrophic model in the cylindrical domain 
O: 



eu- fv+px = 0, (4) 

ev + f u + py = 0, (5) 

Sw +pz ^T, (6) 

Ux + Vy + Wz ^ (7) 

Kh {Txx + Tyy) — KyTzz + uTx + vTy + wTz=Q, (8) 



where e and 6 are positive constants representing the linear (Rayleigh friction) damping coefflcients, and 
Kh is positive constant which stand for the horizontal heat diffusivity. We partition the boundary of SI 
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into: 

T^ = {{x,y,z)€Tl:z = 0}, (9) 

n = {{x,y,z)GTi:z = -h}, (10) 

= {{x, y,z)en: {x, y) e dM, -h<z<0}. (11) 

System (4)-(8) is equipped with the following boundary conditions - with no-normal flow and non-heat 
flux on the side walls and the bottom (see, e.g., [14], [15], [18], [22], [23], [26], [27], [28]): 

dT 

on r„ : = 0, — + aT = 0; (12) 
oz 

dT 

on Tfo : w = 0, ^ = 0; (13) 
oz 

dT 

onT, : {u,v)-n = 0, ^ = 0, (14) 

on 

where ft is the normal vector to Fg. In addition, we supply the system with the initial condition: 

T{x,y,z,0) ^To{x,y,z). (15) 

In this paper we focus on the question of, and prove, the global regularity and well-posedness of the 
3D Salmon's PG model (4)-(8) for all time and all initial data. We remark that a general discussing 
concerning the nonlinear system (4)-(15) was presented in [31], but without providing any evidence of its 
global regularity, a problem that we provide a positive answer for it in this contribution. 

The paper is organized as follows. In section 2, we introduce our notations and recall some well-known 
relevant inequalities. In section 3 we show the short-time existence of strong solutions of system (4)-(8) 
employing a Galerkin approximation procedure. Section 4 is the main section in which we establish the 
required estimates for proving the global existence and uniqueness of the strong solutions, and also show 
their continuous dependence on the initial data. 

2. Preliminaries 

Let us denote by L'^{fl) and W"^''^ (O.) , H'^ {il) the usual L''— Lebesgue and Sobolev spaces, respectively 
(cf., [1]). We denote by 



UWr = { I \(p{x,y,z)\'' dxdydz ] , for every G L'-(fl). (16) 

We set 



~ — dT 

V = TeC~(17): — 



z= — h 



0; i^+c^T 



z=0 



0:% 

on 



and denote by V the closure spaces of V in H^{Q,) under the -ff^ —topology. For convenience, we also 
introduce the following equivalent norm on V: 

Ufv = i^h\\d^4>{x, y, z)\\l + Kh\\dycf>ix, y, z)\\l + {\\d,(l>{x, y, z)\\l + a\\<i){z = 0)||i2(M)) • (17) 

The equivalence of this norm on V to the i?^— norm can be justified thanks to the Poincre inequality (21), 
below. 
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Next, we recall the following three-dimensional Sobolev and Ladyzhenskaya inequalities (see, e.g., [1], 
[9], [10], [13]) 



L3(n) < Co 

Le{Q) < C'o||V'||/i-i(n) 



1/2 

L2(n) 

1/4 



1/2 

3/4 



(18) 

(19) 
(20) 



for every tp G H^{Q). Here Cq is a dimensionless positive constant which might depend on the shape of M 
and O but not on their sizes. We also introduce the following version of Poincare inequality 



(21) 
(22) 



By solving the linear system (4)-(6) we obtain 

^Px + fPy 



u = —- 



w = 



■Pz 



Thanks to (7) we have 



tPx + fPy 

e2 + /2 



+ 



~fPx + epy 



+ 



T 



= 0. 



Using the boundary conditions (12) and (13) we infer the following boundary conditions: 



on r„ and : Pz — T = 0, 



and on Fg : — = 0, 
oe 



(23) 
(24) 
(25) 

(26) 
(27) 



where e = £2i±Zii£i, and k is the unit vector of vertical direction. Notice that by following the techniques 

developed in [12] and [37] (for the case of smooth domains, see, for example, [13] p. 89, and [33]), the 
three-dimensional second order elliptic boundary-value problem (26)~(27) has a unique solution for every 
given T; moreover, this solution enjoys the following regularity properties. Taking the L^{fl) inner product 
of equation (26) with p, integrating by parts and applying the boundary conditions (27) and using the 
Cauchy-Schwarz inequality, we obtain 



■pI) 



Pt 
5 



If 1 
dxdydz = -r / Tp^ dxdydz < t||T'||2||p2||2- 
Jn 



Denote by 



< Fo = min f < Fi = max /. 



(28) 



(29) 



We observe that the assumption Fq > indicates that we are dealing with a mid-latitude case and away 
from the equator. By using (29) and applying Young's inequality to (28), we reach 



/ 

Jn 



{pI+pD + 



Pz 

26 



dxdydz < 



j 

Jn 



e^ + P 



{pI+pD + 



Pz 

26 



dxdydz < 



26' 



12- 



(30) 
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Furthermore, by (26) and the above estimate, we have 



< 



< C 



^/^ /3(lbx||2 + |b,||2) 



/5Px(e'-/') + 2e/3/pj, , T, 



(e2 + /2)2 



TII2 + 



(31) 
(32) 



^eiW2(e2 + ^2)h' ■ J 
As a result of the above and (23)-(25), we obtain 

huh + hvh + \\5wh < CiWVph + \\Th) < C\\Th, 
and 

hu\\Hi(n) + \\ev\\Hi(n) + WSwWhhq) < C{\\Wp\\m{n) + \\T\\m{n)) < C\\T\\H^^l)■ (33) 

Definition 1. Let Tq S V, and let T be a fixed positive time. {u,v,w,p,T) is called a strong solution of 
(4)-(8) on the time interval [0,7^ if 

1) 

TeC{[o,nv)nL\[o,T],H\n)), 

Tt{z = 0) GL\[0,r],H-'/\M)). 

2) {u,v,w,p) satisfies (23)- (27). 

3) Moreover, (8) is satisfied in the weak sense, namely, for every to & [0)T] 



/ T{t)ij; dxdydz - I T{to)tp dxdydz 
Jn Jn 



to 



{KhTx-ipx + KhTyipy + KvT^ipz) dxdydz + k„ a / T{z = 0)'ip{z = 0) dxdy 



M 



ds 



(34) 



/ [ [v- VT(s) + wT^{s)] tp dxdydz ds = [ [ Qip dxdydz ds 
Jtn Jn J to JO, 



for every tp gV, and t G [to,7^- 



3. Short-time Existence of the Strong Solutions 

In this section we will show the short-time existence of the strong solution of system (4)-(8). 

Theorem 2. Let Q £ L'^i^l) and Tq £ V be given. Then there exists a strong solution {u,v,w,p,T) 
of system (4)-(8) on the interval [0,T***], where T*** is a positive time given in (54), below. Further- 
more, dtT e L^{[0,T***];L^{i})) and dtT{z = 0) G L'^{[0,T***]; H-'^/^iM)); and equation (8) holds as a 
functional equation in L^([0,T***]; L^(fi)). 



Proof. We will use a Galerkin like procedure to show the existence of the strong solution for system 
(4)-(8). First, we will show the existence of the weak solutions. Let {^fc e F fl H'^{^)}'^^i and {Afc e 
R+l^-^ be the eigenfunctions and their corresponding eigenvalues of the second order elliptic operators 
—Kh {Txx + Tyy) — KyTzz, subject to the boundary conditions (12)-(14) (see, e.g., [13]). The eigenvalues are 
ordered such that < Ai < A2 < • • • ; moreover, {(^fe}^^ is an orthogonal basis of L^(0). Let m e Z+ be 
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fixed and Hm be the linear space generated by {4>k}^=i- We will denote by Pm ■ — > Hm, the orthogonal 
projection in L^. The Galerkin approximating system of order m that we use for (4)-(8) reads: 

eUm- fVm + da:Pm = 0, (35) 

eVm + fUm + dyPm = 0, (36) 

6w,n+dzPm=Tra, (37) 

dxUm + dyVm + dzWm = (38) 

^t^m (^xx-^m ~t~ ^yy^m) ^v^zzPm ~^ Pm {^m^x^m ~t~ ^m^yPm ~t~ ^ra^zPm\ PmQ ■> (*^^) 

Tm{x,y,z,0) = PmTo{x,y,z), (40) 

where = Y^^=i^k{t)(l>k{x,y, z), and {um,Vm,Wm,pm) is the solution of the system (35)-(38) under 
boundary condition Wml^^o ~ ''^"ilz=-h ~ ^' (^mj^m) ' ^Ir ~ ^- Based on discussion in the previous 
section, equation (39) is an ODE system with the unknown ak(t),k = 1, - ■ ■ ,m. Furthermore, it is easy to 
check that the vector field in equation (39) is locally Lipschitz with respect to ak{t),k = 1, - ■ ■ , m, since it 
is quadratic. Therefore, there is a unique solution ak{t),k = I, - ■ ■ , m, to equation (39) for a short interval 
of time [0,7^]. Let [0,7^*) be the maximal interval of existence for system (35)-(40). We will focus our 
discussion below on the interval [0, 7^*), and will show that 7^* = +oo. 
By taking the L^(f2) inner product of equation (39) with T^, we obtain 

1'3I^ + {WdxTrr^Wl + WdyTrr^Wl) + (||5,T„||i + a||T„(^ = 0)||^) (41) 

+ WmdzTm] Tmdxdydz = / QTmdxdydz. (42) 
It is easy to show by integrating by parts and by using the relevant boundary conditions (12)-(14) that 

[umdxTm + VmdyTm + Wm^^T^] T„ dxdydz = 0. (43) 



/ 

Jn 



In 

Furthermore, by the Cauchy-Schwarz inequality and (21) we have 



/ 

Jn 



QTm dxdydz 



< WQUTmh 



< ^110112 [Kh iWdxTmWl + WdyTmWl) + Kv {\\dzTm\\l + a\\Tm{z = 0)\\l)] , 

V Ai 

where Ai is the first eigenvalue discussed above. From the above estimates, we obtain 

^^^+«h(||a.r„||i + ||5^T„||2)+K„(||a,r„||2+a||T„(z = 0)||i) < (44) 
Consequently, we have, 

^^ + Ai||T„||2<^. 
Thanks to Gronwall inequality, we conclude that 



12 < llToll^ e--^ * + (45) 



for every t G [0,7^")- From the above, we conclude that T^(t) must exist globally, i.e., 7^* = +oo. 
Therefore, for any given T > and any t € [0, 7], we have 



Tmml<\\To\\le-'^'+^^. (46) 
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Furthermore, by integrating (44) with respect the time variable over the interval [0,t], for t £ [0,T], and 
by (46), we get 



y*K(||a,T„(5)||2 + ||5yT„( 
WQWlt 



s)\\l) 





< ll^o||i + 



■ {\\d,Trr,{s)\\l + a\\Trr,{z = 0){s)\\l)] ds 



A? 



(47) 



As a result of all the above we have established that exists globally in time, and that it is uniformly 
bounded, with respect to m, in the L°°{[0,'T]; L'^{fl)) and L'^{[0,'T];V) norms. 

Next, and similar to the theory of 31? Navier-Stokes equations (see, e.g., [9] and [30]), let us show that 
dtTjn is uniformly bounded, with respect to m, in the L^{[0,T];V') norm, where V is the dual space of 
V. From (39), we have, for every tp G V 

[UmdxTm + VmdyTm + WmdzTm] , V) • 

Here, (•, •) is the dual action of V. It is clear that 

|(PmQ,V')l < IIQIUII^Ib, (48) 
and by integration by parts and using boundary condition (12) (14), we have 

\{i^h{d,,Tm + ^yyT„^) + K,d,,Tm,^p}\ < C\\T^\\v llV'lk, (49) 
recall that || • \\v is defined in (17). Next, let us get an estimate for 

\{Pm [UmdxTm + VmdyTm, + WmdzTm] , 1p)\ 

[UmdxTm + VmdyTm + WmO^Tm] tpm dxdydz , 



J 

Jn 



where ipm = -PmV'- Thus, by integration by parts and using (38), (32), (33) and relevant boundary condi- 
tions, we obtain 

I {Pm [UmdxTm + VmdyTm + WmdzTm] , V) I 
= / [Umdxtpm + Vmdytpm + ^m^zV'm] Tm dxdydz 



< C[\\Um\\4. + \\Vm\\i + WWmWi] WTmWi W^i'rnh 



< c 



1/4, 



3/4 



I 1^ 1 1 2 II ''''^'t 1 1 i/i 



1 1/4 1 



J-mll^l T 111^)71112 ll^mllifi T"ll'"m|l2 II ll //i') 11^™ II 2' II II ffi* II ^V'm || 2 
||V^||2. 

Therefore, by the estimates (48)-(51), we have 

mTm,^)\ < C (1IQII2 + ||T„||y + ||T„||2 + ||T„||^/'||T„||^/') 

Thus, we have 

\\d,TmmUt < C (1IQII2 t^/' + \\To\\l + ^Ip) . 



1 1/4 1 



13/4, 



<C\ 



<-m\\2 



\Tm\\l^'\\TnMT 



(50) 
(51) 



(52) 



Therefore, dtTm is uniformly bounded, with respect to to, in the i:3([0,71; V) norm. Thanks to (46), 
(47) and (52), one can apply the Aubin's compactness Theorem (cf., for example, [9], [30]) and ex- 
tract a subsequence {Tm^} of {Tm} and a subsequence {dtTmj} of {dtTm}; which converge to T e 
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L°°{[0,T]; L^{n)) n L^{[0,T];V) and dtT e L^[0,T]:V'), respectively, in the following sense: 

( Tr^.^T in L\[0,n;L^m strongly; 

T„. T in L°°{[0,r];L'^{n)) weak-star; 

T™^. ^ T in L\[0,T];H\n)) weakly; 

t^tT^. ^StT inLt([0,7-];l/') weakly. 

Moreover, from (35) (38) (sec also (4)-(7)) we observe that {urmVm,Wm} depend linearly on Tm- There- 
fore, the elliptic estimates (32) and (33) imply, thanks to (46) and (47), uniform bounds, with respect to 
TO, for {um,Vm,Wm} in ( [0, T] ; (f2)) and L^([0, T]; iJ^(r2)), respectively. Therefore, we can extract 
a subsequence of {umj, VmjWmj}, corresponding to the readily established converging subsequence for the 
temperature {Tm^}, which will be also labeled {umj,VmjWmj}, that converges to {u,v,w} weak-star in 
L°°{[0,T^] L^{Q)), and weakly in L^{[0,7^; H^{Q)). By passing to the limit, one can show as in the case 
of Navier-Stokes equations (sec, for example, [9], [30]) that T also satisfies (34). In other words, T is a 
weak solution of the system (4)-(8). 

By taking the L^(f2) inner product of equation (39) with —Kh{dxxTm + dyyT^) — Kyd^^T^, we obtain 

[Kh {WdxTmWl + WdyTmWl) + Hv {WdzTmWl + a\\T^{z = + \\Kh{d^xTm + dyyTra) + Kyd,,Tm\\l 

= I {Q- UmdxTm + VmdyTm + WmdzTm) (H'hidxxTm + dyyT^) + Kyd^zTm) dxdydz 

< (IIQII2 + WUmhWdxT^h + WVmMdyT^h + WWmMd.T^Ws) \\Khid^xTm + OyyT^) + n^d^.Trah 

< (IIQII2 +C||i;n||6||VT„||3) \\tih{d^xTm + dyyTm) + K^d,,Tmh 

< {\\Qh + C\\Tm\\l) \\Kh{d^xTm + dyyTm) + K^d,,T,nh 

+C [Kh {WdxTmWl + WdyTmWl) + l^v (ll^.T^Hi + a\\Tm{z = 0)\\l)] ^ \\Kh{dxxTm + dyyTm) + tivd,,Tm\\l ■ 

Therefore, applying the Cauchy-Schwarz inequality and Young's inequality to the above estimate, we 
obtain 

^ K (l|a,r„||2 + \\dyT„,\\l) + K„ (||a,r„||^ + a\\T^{z = o)||^)] + |iK,.(a,,T„ + dyyT^) + Kyd,,T^\\l 

< IIQIIi + C\\Tm\\t + C [Kh iWdxTmWl + WdyTrnWl) + {\\d,Tm\\l + a\\T^{z = ml)f ■ (53) 
Consequently, we have 

i^h {WdxTmWl + WdyTmWl) + (||5.T„||i + a||T^(^ = 0)||i) 

^ Kh {WdxToWl + WdyToWl) + i^v {WdzToWl + ct||To(^ = 0)||i) 

~ {l-Ct (llTolli + IIQIIi) [K^(||a,To||2 + ||a,To||2) + /^,(||a,To||i + a||To(^ = 0)||i)])'/'' 

Therefore, for every t e [0,T***], where 



1 

- AC (dlTollt + llOlli) [Kh {WdxToWl + WdyToWD + mT4l + a||To(^ = 0)||i)]) ' 



T*** •= f54l 



we have 



Kh {WdxTmWl + WdyTmWl) + [WdzTraWl + (^\\Tra{z = ^1) 

< 2 [Kh {WdxToWl + WdyToWl) + Kv iWdzToWl + a||To(^ = 0)||^)] . (55) 
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Moreover, by integrating (53) we obtain 

/ \\Kh{dxxTm{s) + dyyTm{s)) + Kydz;,Tm{s)\\l ds 
Jo 

< Kh {\\d..To\\l + lia.Toli^) + (Wd^ToWl + a\\Toiz = 0)\\l) + \\Q\\lt + C (^\\To\\l e"^^ * + t 

+C[K;,(||a,ro|!2 + ||9,To|l2)+«„(||a,ro||2 + a||To(z = 0)||2)]'^ t, te[0,T***]- (56) 

Notice that exists, globally. What we have just proved is that the L'^{[0,T***]; H'^i^l)) norm of 
Tm is bounded uniformly with respect to m. As a result of all the above we have exists, at least, 
on [0,r***] and is uniformly bounded, with respect to m, in L°°([0, T***]; F) and L'^{[0,T***]; H^ifl)) 
norms. Furthermore, and as for the theory of the Navier-Stokes equations (see, for example, [9], [30]), 
we can use the above bounds (55) and (56) to show that the L^([0, T***]; i^(0)) norm of dtTm and the 
-L^([0,T***];i?~"^/^(M)) norm of dtTjn{z = 0) are uniformly bounded with respect to m. Passing to the 
limits, we conclude that there is a strong solution to system (4)-(8), at least, on [0,T***]. Furthermore, 
this strong solution enjoys the following properties: 

9tr e l2([o,7-***];L^(0)) and dtT{z = 0) e L^{[{),r***]:H-^/'^{M)). (57) 

The above regularity estimates are sufficient to complete the proof of Theorem 2, following standard 
techniques from the theory of the Navier-Stokes equations (see, e.g., [9] and [30]). Furthermore, as a 
consequence of the above estimates, in particular those implying (57), we conclude that equation (8) holds 
as a functional equation in ^^([0, T***]; L^(f^)). □ 

4. Global Existence and Uniqueness of the Strong Solutions 

In the previous section we have established the short-time existence of the strong solution to system 
(4)-(8). In this section we will show the global existence and uniqueness, i.e. global regularity, of strong 
solutions to the system (4)-(8), and their continuous dependence on initial data. 

Theorem 3. Let Q G L'^{^), Tq G V and T > 0, be given. Then there exists a unique strong solution 
{u,v,w,p,T) of the system (4) (8), on the interval [0,7^, which depends continuously on the initial data 
in the sense specified in equation (72) below. 

Proof. Denote by {u,v,w,p,T) the strong solution corresponding to the initial data Tg with maximal 
interval of existence [0,7^), that has been established in Theorem 2. We will show that = oo. To show 
this we assume by contradiction that % < oo. Consequently, it is clear that 

limsup \\T{t)\\m(n) = oo, 

becausc, otherwise, and by virtue of Theorem 2, the solution can be extended beyond the maximal time of 
existence, T*. Next, we will show that ||T(i)||fl'i(Q) is bounded uniformly on the interval [0,%). In what 
follows we will focus our discussion and estimates on the finite maximal interval of existence [0,%). 

4.1. estimates. As a result of Theorem 2, equation (8) holds in -^jQp([0, 7^); L'^i^)), therefore we can 
take the inner product of equation (8) with T, in L^(il), and obtain 

+ {\\dxT\\l + \\dyT\\l) + (lia.rll^ + a\\T{z = 0)||i) 

= / QTdxdydz- / {ud^T + vdyT + wd^T)T dxdydz. 
Jn Jn 
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After integrating by parts we get 

{uda^T + vdyT + wd^T) T dxdydz = 0. (58) 



/ 



m 

As a result of the above we conclude 

+ (II^.Tlli + ||a,T||i) + {\\d,T\\l + a\\T{z = 0)||i) 

= / QT dxdydz < WQh \\T\\2. 
Jn 

Using (21) and the Cauchy-Schwarz inequality we obtain 

d\\T"^ 



dt 



2kh {\\d.T\\l + WdyTf^ + (||5.T||^ + a\\T{z = (59) 



<^- + -j||Q||^. (60) 
By the inequality (21) and thanks to Gronwall inequality the above gives 

||T||2 < e" ||To||2 + (2/^2 + 2h/af ml, (61) 

for are t G [0,T*). Moreover, we also have 

r K (11^-1^112 + WdyTWl) + {\\d,T\\l + a\\T{z = 0)\\l)] ds 
Jo 

< 2{h^ + -)\\Q\\l t + e- ^TOA^I IITolli + (2/1^ + 2h/af\\Q\\l (62) 



a 



for are t G [0,%). 

We remark that estimates (61) and (62) also follow directly from (46) and (47), respectively. 

4.2. estimates. Recall from Theorem 2 that T e L^^{[0,%), H^{n))r]Lf^^{[0,%), H'^in)), therefore 
|T|^T e L^^^([0,7;);i^(f2)). Since by Theorem 2 equation (8) holds in L^^^{[0,%); L'^{0.)) we can take 
the inner product of the equation (8), in 1/^(0), with |T|*T to get 

+ {Wd.nl + \\dyT\\l) + K4d,T\\l] |T|4 dxdydz + a^4T{z = 0)\\l 

= / Q\TfT dxdydz- / {uT^ + vTy + wT;,)\T\^T dxdydz. 
Jn Jn 

By integration by parts, and using (7) and the boundary conditions (12)-(14) we get 

{uT^ + vTy + wT^) \T\'^T dxdydz = 0. (63) 



/ 

Jn 



In 

As a result of the above we conclude 



I'^ll^i+S / [Kh{\\d,T\\l + \\dyT\\l)+K4d.T\\l] \T\Uxdydz + aK4T{z^0)\\l 



6 dt 



= [ Q\T\'T dxdydz < WQhWnlo < C\\Qh {\\T\\l HVT^H + \\T\\l) . 
Jn 
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By the Cauchy-Schwarz inequality we get 



dt 



16^ 



/ [nh{\\d,T\\l + \\dyT\\l) + K^d.Tf^ \T\Uxdydz + aK4T{z = m\ 
Jn 



= / Q\T\^T dxdydz < CWQWiWnt + \\Qh\\T\\l < C||g||^||T||| + ||T||t 
Jq 

Thus, from the above and (22), we have 

^ < CWQWl + \\T\\l < C [ml + \\T\\l + [Wd^Ml + WdyTf^ + '^.II^.TII^] . 
By integrating the above inequahty and using (61) and (62), we get 

||T(t)||^ < C [(1 + IIQIli) {l + t) + ||To||l,.(o)] . (64) 

4.3. estimates. Recall again that T e L^^^{[Q,T*),H^{n))r\ L'^^^{[Q,%),H'^{^)), and since, by The- 
orem 2, equation (8) holds in Lj^^([0, Tli,); ^^(^1)) we can take the inner product of the equation (8) with 
—Kh {Txx + Tyy) — KyTzz, iH ^^(fJ), BxiA use (57) to obtain, thanks to a Lemma of Lions-Magenes concerning 
the derivative of functions with values in Banach space (cf. Chap. III-p.l69- [30]), 

~ K (11^-1^112 + WdyTWl) + (||5.T||^ + a\\T{z = 0)||^)] + \\nn (T.. + T,,) + n^T^^l 

= - / Q{Kh (Txx + Tyy) + KvTzz] dxdydz + / {ud^T + vdyT + wd,T) [kh (T^x + Tyy) + k„T^^] dxdydz 

< [\\Qh + (Me + IHe + \\wh) IIVTIIs] \\Kh (Txx + Tyy) + KyT,,\\^ 

< \\\Qh + CWTWl^^ \\Kh (Txx + Tyy) + «,T,,||^/'" 

\\l^h {Txx "I" Tyy) + KljTj2:||2 • 

By the Cauchy-Schwarz and Young's inequalities we obtain 

^ [Kh i^dxTWl + \\dyT\\l) + {\\d,T\\l + a\\T{z = 0)||^)] + \\Kh {Txx + Tyy) + k,T,,\\1 
<C\\Q\\l + C\\T\\l 
By Gronwall, we get 



{\\d.T{t)\\l + \\dyT{t)\\l) + {\\d.T{t)\\l + a\\T{z = 0){t)\\l) 

+ I \Wh{Txx{s) + Tyy{s)) + K„T,,{s)\\l ds 
Jo 

<C{l + \\Q\\l + \\Tfe}t+\\To\\HHn) 



3 



Thus, 



< C(l + ligil^) t + C {1 + WQWi) {l+t) + llToll^i(^) t + IIToIIhmo) =: K,{t). (65) 



limsup||T||jfi(n) =/s:^(r,). 



This contradicts the assumption that % is finite, therefore, % = oo, and the solution {u,v,w,p,T) exists 
globally in time. 
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4.4. Uniqueness of the strong solution and continuous dependence on initial data. Next, we 
show the continuous dependence on the initial data and the the uniqueness of the strong solutions. Let 
{ui,Vi,Wi,pi,Ti) and (^2, U2, W2,P2, ^2) be two strong solutions of the system (4)-(8) with corresponding 
initial data (ro)i and (ro)2, respectively. Denote hy u = ui — U2,v = vi — V2,w = wi — 'W2,p = Pi — P2 
and e = Ti - T2. It is clear that 

eu- fv+px = 0, (66) 
€v + fu+py = 0, (67) 
6w+p,^9, (68) 
Ux + Vy + Wz = (69) 
dtO - Kh (6i^x + Oyy) - i^vOzz + Ui6'^ + v^Oy + Wi6z + ud^^Ti + vdyT2 + wdzT2 = , (70) 
and {u,v,w) and 6 satisfy boundary conditions (12)-(14). By Theorem 2 and Theorem 3 equation (70) 
holds in L'^{[0,T];L'^{n)) and e L'^{[0,T), H^{n)) n L'^{[0,T), H'^{n)), for all T > 0. Therefore, by 
taking the inner product of equation (70) with in -L^(f2), and using boundary conditions (12)-(14), we 
get 

ld\\0\\? 



2 dt 



{WdMl + \\dyO\\l) + K^WdMl + aPiz = 0)111 



= - / [uiO^ + vi0y + wiOz + u{T2)x + v{T2)y + w{T2)z] dxdydz. 

By integration by parts and again boundary conditions (12)-(14), we get 

- / [ui6x + vi0y + wiOz] dxdydz = 0. (71) 
Jn 



Notice that 



/ [u{T2)x + v{T2)y + w{T2)z] dxdydz 

JQ 

< c\\VT2\\2 (hWl^'Wuf^t + Ml^'Wvf^t + Ml^'Wwf^t) \\e\ 

< c\\^T2\\2\\e\\l^'\\e\f^? < CIIVT2II2 (ii^iii + ml^'iMwl^') . 



<C||VT2l|2(h||4+|k||4+||w||4) 



Thus, 



+ {\\dM + \\dyO\\l) + K^„\\ml + c^Mz = 0)11^ 



<c!|vr2||2(||^||^ + ||0||^/'||V0 

By Young's inequality, wo get 

d\\0r 



If 



d.OWl + \\dyO\\l) + K4dz9\\l + apiz = 0)\\l 



dt 

<c\\vT2\m\i 

Thanks to Gronwall inequality, we obtain 

||^^(^)||2< ||^?(t = o)i|2e^/oiivT.wii|ds^ 

Since T2 is a strong solution, we have by virtue of (65) 

\\0ml<\m = O)\\le''foK%s)ds^ (72) 

where the value of Tq in the definition of Ky in (65) is replaced by T2(0). As a result, the above inequality 
proves the continuous dependence of the solutions on the initial data. In particular, when 0{t = 0) = 0, we 
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have 9{t) = 0, and consequently also u{t) = v(t) = w(t) = 0, for all t>Q. Therefore, the strong solution 
is unique. 

□ 
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